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Both nonmagnetic and magnetic impurity effects in spin singlet and triplet excitonic insulators 
were investigated. The bound state energies caused by single impurity were given. The different 
compositions of the bound states can be used to detect the symmetry of the excitonic insulators. 
In finite concentration problems, nonmagnetic impurities showed same pair-breaking effect in sin- 
glet and triplet excitonic insulators while magnetic impurities showed weaker pair-breaking effect in 
triplet excitonic insulators than in singlet ones. The pair-breaking effects suppressed the ferromag- 
netic range via doping and gave a natural explaination for experimental results. 
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I. INTRODUCTION 

The excitonic phase has been researched for a long 
time since I960—"—. It was proposed in semimetals, in 
which there is a band overlap and the normal states are 
unstable. The true ground state becomes a state with 
electron-hole pairs due to the Coulomb interaction. This 
kind of state can also occur in semiconductors once the 
energy gap becomes smaller than the binding energy of 
an excitoni. This so-called 'excitonic insulator' has been 
discussed for many years theoretically but materials in 
experiment were seldom reported so far. IT — TiSe2, 
among the transition-metal dichalcogenides, may be one 
example of the excitonic insulators which obtained more 
and more direct evidences from experiments^!^. Another 
possible example is CaBg. The latter one has gained 
more attention since weak ferromagnetism was found in 
shghtly doped CaBe in 1999^. Theoretically many au- 
thors relate this to the historical excitonic insulators ^"^. 
By now there are still debates on this proposal, espe- 
cially about the band gap. The angle-resolved pho- 
toemission spectrum showed a gap of le v^^i^^ which 
is too large for excitonic instability and another large 
gap of 0.8ev is given by the LDA+GW calculations^^. 
Band overlap and semimetal framework were demon- 
strated by the de Haas-van Alphen and Shubnikov-de 
Haas experiment&i2r,i§. and local density approximation 
band structure calculations-!^-"-^. Though there are con- 
troversies on the band structure, the excitonic insulator 
scenario is nevertheless a very promising explaination for 
the main physics in CaBg. 

In Cai_a;La2:B6, slightly La doping would induce not 
only the carriers but also impurities / disorders inevitably 
and the impurity effects played an important role in the 
formation of the weak ferromagnetism. According to the 
first principle calculations, the defects or vacancies in the 
boron sublattice can result in impurity bands or local- 
ized gap states and this modification may be important 
to the origin of ferromagnetism in CaBgr^^. From the 
language of excitonic insulators, the impurities will in- 
duce in-gap bound states and it was confirmed by tunnel- 
ing measurements on high-quality defect-controlled single 



crystals2i. Therefore the impurity effect can not be ne- 
glected in the excitonic insulator theory for CaBg. To 
explore how the charged impurity doping can affect the 
excitonic insulators constitutes one of the motivations 
of this paper. Another motivation originates from the 
incompleteness of the research on impurity effects in ex- 
citonic insulators. Historically the impurity scattering 
in excitonic insulators was investigated by J. Zittartz in 
1967— and in that work only intra-band nonmagnetic 
impurity in exciton insulators with no regard for the 
spin freedom was studied. Inter-band nonmagnetic im- 
purity scattering was investigated by T.Ichinomiyai^ and 
other works concentrated in the bilayer system a^^^i^^ in 
which the spin freedom was also neglected. In this paper 
we will discuss and compare the inter-/intra-band mag- 
netic/nonmagnetic impurity effect in singlet and triplet 
excitonic insulators. We will also modify the exciton in- 
sulator scenario on ferromagnetism in Cai_a;Laa;Bg. 

This paper is constructed as follows: In the first sec- 
tion we will introduce our model and Hamiltonian. Then 
we will give the locations of the bound states induced by 
single magnetic or nonmagnetic impurity. Finite concen- 
tration problems of impurities and pair breaking effects 
will be discussed in the third section. In the fourth sec- 
tion we will investigate the impurity caused-modification 
to the ferromagnetism in doped exciton insulators. We 
will leave the summary and further discussion to the last 
section. 



II. MODEL AND HAMILTONIAN 

We introduce the noninteracting Hamiltonian in a two 
band system: 

here the kinetic energy in two electron bands are e^. = 
1^ - \Eg, Efc = + hEa, in which Eg is the band 

gap or overlap and \Eg \ should be quite small to keep the 
ground state remain in the excitonic phase. Without los- 
ing generality, we consider two symmetric bands, which 
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means rua = mi, = m. Now we introduce the interaction 
term, which reads: 



Hi 



Ui 



E 

k.k' ,q,a,(T' 



k+q,cr^ 



where Ui represents the interaction term which comes 
from the Coulomb interaction. In most of the bulk sam- 
ples, the system will choose singlet or triplet excitonic 
phase as the ground state^^ because there is a finite split- 
ting between singlet and triplet states. The splitting may 
come from short-range Coulomb terms which favor triplet 
state, or electron-phonon interaction which favors singlet 
state^^. In the presence of nesting effect, if the screened 
Coulomb interaction term dominates there will be a de- 
generacy between singlet and triplet excitonic phase. In 
this case there is an SU(2)xSU(2) symmetry in the ex- 
citonic insulators. In bilayer systems this splitting is ex- 
ponentially small because of the inter-layer distance and 
the barriers between two layers^*, and spin freedom can 
thus be neglected. It is also believed that the splitting in 
CaBg is quite smallS.. 

Now we introduce the singlet and triplet order param- 
eters: As and At and rewrite the mean-field Model 
into Nambu space: By introducing the vector ip — 
(cLk-fj 0-ki: bk^i bki), the Hamiltonian can be written as : 



Hmf = J2 ^^^^^ 



where M is a 4 x 4 matrix, then the bare Green function 
for singlet excitonic insulators is 



G^\k,uj)=uj- 
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and triplet excitonic insulators: 
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The impurity scattering term can be written as: 

Jin ■ I + Jlm<^ ■ S J2n + J2m<^ ' S 



V 



J2n + J2m<^ ' S Jin ' / + Jlm<^ ' S 



in which h' — 71,111, represents the nonmagnetic and mag- 
netic impurity scattering. S denotes the spin operator of 
the magnetic impurity and z = 1,2 represents the intra- 
band and inter-band channel, respectively. 



III. SINGLE IMPURITY PROBLEMS 

To obtain the bound state energy, one has to get the 
T-matrix which can be written as: 



T 
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1 - Go(r = 0,u)V 



the pole of the T-matrix gives the bound state energy. In 
singlet excitonic insulators, for nonmagnetic impurity. 
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where ai„ = JinirNp , a2r. 
netic impurity case. 
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A5[l - (qi„, + a2mf] 

1 + (aim + a2m)^ 

-As[l - {airn - a2r,if 



1 + (aim - a2m)^ 

where aim = Jim'^NFS{S+l) , a2m = J2m'^NFS{S+l). 
Notice some extreme cases in the singlet excitonic insula- 
tors, for the negative order parameter and if there is only 
one scattering channel exists, the bound state energy will 
have the same form of the single magnetic impurity in s- 

wave superconductors^^— , which reads Eb ^ 

(i=ln,2n,lm,2m). In this case the bound state energy 
will move to the Fermi energy with the increasing of the 
scattering constant while move to the gap-edge side with 
the decreasing of impurity scattering. If there is no im- 
purity scattering (J^^ = 0), the excitation will emerge 
into the normal excitonic gap Eb = ±A. 

In triplet case, the bound state energy can be obtained 
in the same way. For single nonmagnetic impurity 

At[1 - (ai„ + a2n)^] ^A /n 

Ebi = ± r~TZ ~~9. ; ±(ain + a2n)AT < 



Eb2=± 



1 + (ai„ + a^„ 
At[1 - (ai« - a2n)^] 
1 + (ai„ - a2,i)^ 



±(ain - a2n)AT < 



Since the form of the bound state energy for single mag- 
netic impurity in triplet excitonic insulators is quite com- 
plex, we only give the numerical results. The bound state 
energies for single nonmagnetic and magnetic impurity in 
singlet and triplet excitonic insulators are shown in Figl. 

From Figl we can see that location of the bound state 
energy which is induced by single impurity is very differ- 
ent between singlet excitonic insulators and triplet ones. 
Furthermore, in excitonic insulators, the quasiparticle in- 
duced by single nonmagnetic or magnetic impurity can 
consist one or two hole (w < 0) or one or two particle 
type (a; > 0) states in some range of Ji„, J2„ or Jim, 
J2m- This is quite different from the single magnetic 
impurity in s-wave superconductors in which hole and 
electron type excitations must appear symmetrically at 
the same time. It is well known that both the hole- and 
electron-type components of the bound states can be de- 
tected by the low-temperature tunneling microscope^. 
Thus the results of the bound state energies are valuable 
since they can be used to detect the symmetry of the ex- 
citonic insulators, especially for the materials with qutie 
small splitting of singlet and triplet states. 
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Figure 1. Bound state energy induced by single impurity 
in singlet and triplet excitonic insulators via the scatter- 
ing strength. From top to below is the case for nonmag- 
netic/magnetic impurity in singlet excitonic insulators, non- 
magnetic/magnetic impurity in triplet excitonic insulators, 
respectively. Please note that in the figure solid symbols rep- 
resent As,T > and open symbols represent As,t < 0. 



Another interesting finding is that magnetic impurity 
in triplet excitonic insulator is quite different from other 
cases. We can see the bound states caused by single mag- 
netic impurity in triplet exciton insulators are very close 
to the gap edge. With the increasing of the impurity con- 
centration an impurity band is expanding and the states 
are stretching to the Fermi level. Under this condition 
we could expect that finite concentration of impurities 
should behave a weaker pair breaking effect than in the 
other cases and we shall justify this conclusion below. 



IV. FINITE CONCENTRATION PROBLEMS 
AND PAIR-BREAKING EFFECTS 

To investigate the finite concentration problems we in- 
troduce the renormalized Green Function. For spin sin- 
glet excitonic insulators: 



-As-I i^ + i^k 



1,2 E, 



2m 2 



(3) 

with the full Born approximation one can get the self- 
energy : 



< VG{k,iuj)V >i 



where <>/ means averaging all the impurity positions 



and 



represents the density of impurities. Then we 



can calculate the self-energy which is given by 
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As we know the renormalized Green Function can be 
written as: 

G-\k,iu) = G,\k,u;)-J2i^) 
then we have the self-consistent equations: 



As - 



(4) 
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Note here As represents the excitonic order parameter of 
singlet excitonic insulators in the presence of impurities. 
We consider the one-channel (inter-band or intra-band) 
impurity scattering case, then the impurity effect only 
contributes to the renormalization of frequency and the 
order parameter. Making a; — > ioj, the order parameter in 
zero temperature and the excionic gap can be got from: 



A = VNf 



n = A(i 
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where ec is the cutoff and C = [/(w) — /(A)]/A. In the 
expression of /(w) and /(A) can be well defined as the 
factor of the renormalized frequency and order parameter 
if only one impurity scattering channel opens ( Ji^ = or 
J2^ = 0). They read 



/(A) = 



Jlu or J2(y— 



1 



Ji„ or J2„=0 



The same analysis can be applied to the triplet case, 
in which the renormalized Green- Function can be written 
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The self-consistent equations for one channel impurity 
scattering are: 

,1 1 , w 
w = w + ( , )- 



At 



At = AtT ( , ) 

Ttln Tt2m 



Note here the {+,m) represents the magnetic case and 
(— ,n) represents the nonmagnetic case and At is the 
excitonic order parameter of triplet excitonic insulators 
in the presence of impurities. Then the scattering times 
read: 
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1 
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1 
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The numerical results for the order parameters and ex- 
citation gaps with the impurity doping in singlet and 
triplet excitonic insulators are shown in Fig2. 
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Figure 2. The order parameters and the excitonic gaps via the 
impurity concentration. The left figure is the nonmagnetic 
impurities in singlet and triplet excitonic insulators and the 
right figure is for magnetic impurities case. Aso and Ato 
represent the order parameter of singlet and triplet excitonic 
insulators in the absence of impurities. 

From Fig2 we can see that both of magnetic and non- 
magnetic impurities have the pair breaking effects in ex- 
citonic insulators and there always exists a gapless region 
between the gaped excitonic insulator phase and the nor- 
mal phase. This phase transition from gaped phase to 
gapless phase is of second order as pointed out by J. 
Zittarta^^. However, the nonmagnetic impurities have 
the same pair breaking effect to singlet and triplet ex- 
citonic insulators but magnetic impurities behave differ- 
ently in these two systems. A significant difference can be 




seen from the right part in Fig2 in which the triplet exci- 
tonic insulators are more robust with magnetic impurity 
doping. Notice that the pair breaking effect coming from 
the self-consistent equations is decided by the relative 
amplitude and sign between the factor of renormalized 
frequency and order parameter in the self-energy. There- 
fore, to compare the pair-breaking effect in quantization 
between different cases, we define the pair breaking factor 
which can be written as 

It is analogical to the notations in R. Balian et al 
's work in p-wave superconductors'^'^. According to 
Abrikosov-Gor'kov theory^, in conventional s-wave su- 
perconductors, X = —I stands for magnetic impurities 
and a; = 1 stands for nonmagnetic impurities. Then we 
list all the pair breaking effects in singlet and triplet ex- 
citonic insulators. 



(8) 



where Xs,n, Xs,m, At.m, At,™ represent the nonmagnetic 
impurity scattering and magnetic impurity scattering in 
the intra-band/inter band channel in singlet and triplet 
excitonic insulators, respectively. 

It is well known that the nonmagnetic impurities are 
not pair breakers in s-wave conventional superconduc- 
tors according to Anderson's theorem. Note in singlet 
and triplet excitonic insulators nonmagnetic impurities 
have the same pair-breaking effect. This is analogical to 
the magnetic impurities in the s-wave superconductors in 
which the time reversal symmetry is broken due to the 
spin-flip impurity scattering. The Kramer's degeneracy is 
replaced by the particle-hole symmetry in the singlet and 
triplet excitonic insulators. According to our calculations 
triplet excitonic insulators are more stable with magnetic 
impurity doping than singlet samples, ht/^s — 3, where 
riT and ns are the critical doping concentration of mag- 
netic impurities in triplet and singlet excitonic insulators, 
respectively. This result is in accordance with the single 
impurity analysis. From Figl we can see that the bound 
states energy is close to the gap edge for the magnetic im- 
purity in triplet excitonic insulators. And it needs larger 
impurity concentration in triplet than in singlet excitonic 
insulators to fill the gap with states. 



V. IMPURITY EFFECTS AND 
FERROMAGNETISM 

In undoped CaBg the splitting of singlet and triplet 
excitonic phases is neglectable^, then the SU{2) x SU{2) 
symmetry is preserved and in this case the samples show 
no magnetic signals since there is no polarization in the 
systems. All kinds of spin polarization states have the 
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same energy. Once the doping induce extra carriers 
the degeneracy between singlet and triplet states can be 
lifted^i^. The doped electrons are aligned in one direc- 
tion to make the energy minimized. Then the system 
have two order parameters A-|- and A^, which have the 
relationship with A5 and A^ in doped situation^ 



A;Zn(-^)5 + At/n(AsoATo)^ = A^ln\A^\ 

Ato 



At^"(^)^ + AXAsoAto)^ 



A^ln{2n + J An? + A|) 



(9) 



where A50 and A^o are the singlet and triplet exci- 
tonic order parameters of undoped case in the absence 
of impurities, respectively, n is the concentration of the 
doping electrons. At zero doping at T = we have 
A50 = Ato = Aq and in doped samples we hav o^'^^ 



A^ = v/Ao(Ao 
At = Ao 



An) 



(0 < n < A/4) 
(0 < n < A/2) 



Now we can see that ferromagnetism can exist in the dop- 
ing range of [0, Ao/2]. In the doping level [0, Ao/4], the 
doped electrons are polarized (DP). In the doping level 
[Ao/4, Ao/2], all of the electrons and holes are paired in 
one direction and system is in the fully polarized phase 
(FP). When n > Ao/2, the excitonic state is no longer 
favored and system undergoes a first order transition into 
the normal metallic state. This phase transition is differ- 
ent with the excitonic insulator-metal phase transition we 
discussed above. The former one is caused by carrier dop- 
ing and becomes a normal metal after the transition. The 
latter one is caused by the impurity doping and under- 
goes a transition into the impure metal phase. Besides, 
for the impurity induced excitonic insulator-metal phase 
transition, generally there is a second order phase transi- 
tion into a gapless region before the excitonic insulator- 
metal phase transition. 

In bulk materials the La doping would induce impu- 
rities or disorders inevitably and in this case both the 
carrier doing and impurity effect must be considered. 
For CaBg, since there is no exact experimental data of 
Ao and if we consider a general excitonic gap of order 
^ GOO-fC— , the ferromagnetism range via doping should 
be — 3% and this doping range is much larger than the 
experimental results 0—1%. If the excitonic insulators 
scenario is responsible for the physics in CaBg, one pos- 
sible reason for the narrow range of ferromagnetism is 
the anisotropy of the band structure^'^. However this is 
not enough and we must involve the pair breaking effects 
caused by impurity or disorder. A direct evidence for 
this comes from the tunnelling experiments^^, in which 
it is surprising to find that there exists ferromagnetism in 
99.9% pure boron of CaBg within a small range of carrier 
doping while there are no magnetic signals in 99.9999% 
samples. For high purity samples the exciton gap may be 
too large compared to the doping concentration and the 



splitting between A^ and A^ is neglectable. Therefore 
the ferromagnetism can not be found. More disorders 
or impurities in the samples suppress the excitonic order 
parameters and make it comparable to the doping level. 
Then the impurity effect results in two changes: a more 
significant splitting of A^ and A^ and a narrower doping 
range of ferromagnetism. 

To express our physics in detail we carry out the nu- 
merical work by replacing the A in (1) with As-I+Ax-cr. 
For simplification we choose the triplet exciton along the 
z direction. Then with the same full Born approximation 
we have the self-consistent equations with nonmagnetic 
impurity scattering: 



Wi = a; + 



U)2 — 1^ + 



UJi- 



{As + AtY - Cjf 
{At - As)2 - (i| 



At + A5 = (Ai 



At - As = (Ai 



(10) 

the self-consistent equations with magnetic impurity 
scattering: 



(At + As)^ 

(At + As)2-^2 

(At-As)^ 
(At - As)2 - c^l 
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UJ2 = + 
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(At + As)2 -uf J {At - As)^ - Co. 



UJl- 



At + As = (At + As) 



{At + AsY-CoI ^J {At ~ Asf - CjI 
(AT + As)-i- 



(At + As 
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(At- As): 



(At - As)2 - cif 



(At + As)^ 

At - As = (At - As) + • 



(At + As 



1 2 _ ,7,2 



(Ar-As)^ 
(At - As)2 - 



(11) 

Taking both the impurity effects and carrier doping into 
considerations, we have to solve the Equs.(6), (9) and 

(10) for nonmagnetic impurities and Equs. (6), (9) and 

(11) for magnetic impurities. Now the ferromagnetic 
range and the fully polarized range should be — ni and 
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712 — ni, where the critical charged impurity concentra- 
tion Hi and n2 satisfy 



A- 



imp 
C 

imp 



ni. 



n2) 



Note here the impurity concentration Uimp and density of 
doping electrons should be replaced with charged impu- 
rity concentration nfj^p. Our numerical results are shown 
in Fig3. 
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Figure 3. Order parameters with doping concentration. 
Note in the figure and Uimp represent the charged 

and non-charged impurities concentration, respectively. From 
top to below: _ = 0.5; _ ^ „ , _ ^ = 1 and 



1 1 = 1.5. 



In Fig3 we note that in the nonmagnetic impurity scat- 
tering case, the pair breaking effect is the same for sin- 
glet and triplet excitons. But in the magnetic impurity 
scattering case, the order parameter of singlet and triplet 
excitons splits with the increasing of the impurity concen- 
tration. Both two kinds of order parameters disappear 
at the same doping concentration at which the excitonic 
insulator becomes an impure metal. In the situation here 
we do not consider the gapless excitonic insulator since 
this regain is quite small(< 0.05Ao). From top to be- 
low in Fig3 we can also see that the ferromagnetic region 
— ni is growing narrower with the increasing of the 
impurity scattering. To express that clearer we plot the 
phase diagram in J— n plane in Fig4. It is very clear that 
the nonmagnetic and magnetic impurities would suppress 
the ferromagnetic range via doping. In magnetic impu- 
rity case, the fully polarized range is destroyed by the 



Figure 4. Phase diagram in the doping concentration- 
impurity scattering plane. In the figure IM represent the im- 
pure metal phase. FP and DP represent the fully polarized 
excitonic insulator phase for all the electron-hole pairs and 
all doped carriers, respectively. Note the weak ferromagnetic 
region includes FP and DP region. 



impurity scattering in all doping levels. For pure car- 
riers doping the ferromagnetic range should be — 3% 
and this region are suppressed to — 1% with the im- 
purity scattering JinirNp ~ 2.8 for nonmagnetic case 
and JimTiNp — 2.0 for magnetic case, where i = 1,2. 
And this result is in accordance with the experimental 
data. It is worthy to point out that here we consid- 
ered the charged impurity case which means the density 
of doped electrons and impurities have the same value 
n — Uimp — nf^p. In CaBg, the charged impurities are 
induced by La doping. If there are additional impuri- 
ties/disorders which are not accompanied with electrons 
doping, which means riimp > n, they can strengthen the 
pair-breaking effects and these extra impurity scattering 



terms can be absorbed into the factor of Ji, 
rectly. 



or Jim di- 



VI. DISCUSSION AND SUMMARY 

We investigated the nonmagnetic and magnetic impu- 
rity effects in singlet and triplet excitonic insulators in 
the inter-band and intra-band channels. In single impu- 
rity problems, the bound states showed different compo- 
sitions in singlet and triplet excitonic insulators, which 
can be used to detect the symmetry of excitonic insula- 
tors in experiments. In finite concentration problems, all 
kinds of impurities showed pair-breaking effects. In the 
presence of impurities, the excitonic insulator scenario 
for ferromagnetism in CaBg was modified . The doped 
excitonic insulators can be characterized in three classes 
by different kinds of doping: 

(1) Electrons doping but no impurity effects induced: 
Ferromagnetism emerges with the electrons doping and 
there are two ferromagnetic phases with doping, the po- 
larized phase for the doped electrons and the full polar- 
ized phase for all the excitonic pairs. 

(2) Impurity effects induced but no electrons brought 
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in: i. Nonmagnetic impurities in one channel show same 
pair breaking effects to singlet and triplet excitonic in- 
sulators, ii. Magnetic impurities in one channel show 
weaker pair breaking effect to triplet excitonic insulators 
than singlet ones. 

(3) Both electrons and impurity effect were induced 
by doping: Nonmagnetic and magnetic impurity effects 
suppressed the ferromagnetic region and this provides an 
natural explaination for the quite narrow region for fer- 
romagnetism in the experiments of CaBg. In the mag- 
netic impurity case, the fully polarized phase for all the 
electron- hole pairs was destroyed in all doping levels. 

The model we discussed here is based on two nested 
bands which have the same effective masses. Once apart 



from this situation, the two bands have different effec- 
tive masses, nia ^ ruh- In this case the impurities are 
still pair breakers because they can be seen as oppositely 
charged pairs to an excitoni^^ . This may make the quan- 
titative calculations much more difficult but should show 
the same conclusions as ours qualitatively. The more 
anisotropy of band structure will obviously reduce the 
total magnetism and may lead to the well-known Larkin- 
Ovchinnikov-Fulde-Ferrerll state^i^ which is not consid- 
ered in this work. 
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